Next-to-Leading Order Hard Scattering Using Fully Unintegrated Parton Distribution 
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We calculate the next-to-leading order fully unintegrated hard scattering coefficient for unpo- 
larized gluon-induced deep inelastic scattering using the logical framework of parton correlation 
functions developed in previous work. In our approach, exact four-momentum conservation is 
maintained throughout the calculation. Hence, all non-perturbative functions, like parton distri- 
bution functions, depend on all components of parton four-momentum. In contrast to the usual 
coUinear factorization approach where the hard scattering coefficient involves generalized functions 
(such as Dirac (5-functions), the fully unintegrated hard scattering coefficient is an ordinary func- 
tion. Gluon-induced deep inelastic scattering provides a simple illustration of the application of the 
fully unintegrated factorization formalism with a non-trivial hard scattering coefficient, applied to 
a phenomenologically interesting case. Furthermore, the gluon-induced process allows for a param- 
eterization of the fully unintegrated gluon distribution function. 

Keywords: QCD, factorization 



I. INTRODUCTION 



The standard factorization theorems of perturbative 
QCD (pQCD) are the main ingredients in many phe- 
nomenological calculations of high energy processes (for 
a review of the usual approach to colhnear factorization, 
see [l|). They are especially important in Monte Carlo 
event generator (MCEG) calculations (see e.g., 0]) and, 
more generally, in many of the calculations relevant to 
up-coming experiments such as those taking place at the 
Large Hadron Collider. The standard derivations, how- 
ever, rely on a number of approximations that alter final 
state momentum such that over all four-momentum con- 
servation is violated. Hence, there is a potential for large 
errors when the standard factorization theorems are ex- 
tended to more general situations than what have been 
considered in the past, particularly when the details of 
final states are of interest. 

To understand the origin of these issues, it is helpful 
to begin by briefly reviewing the basic structure of the 
more standard approaches to factorization. To start off, 
we recall that a basic component of the standard colhnear 
factorization theorems is the concept of a universal par- 
ton distribution function (PDF) which describes the like- 
lihood of finding a parton in the proton with a particular 
plus component of momentum. The partons are treated 
as having zero transverse and minus components of mo- 
mentum in the hard scattering calculation. In the stan- 
dard colhnear factorization theorems, the PDFs have pre- 
cise definitions as expectation values of appropriate com- 
binations of field operators [S, 3- In these definitions, 
there is an integration over the small transverse and mi- 
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nus components of parton momentum because they are 
neglected in the hard scattering calculation. Thus, the 
standard PDF may be referred to as the "integrated" 
PDF. Standard colhnear factorization relies directly on 
the approximations that allow the small components of 
four-momentum to be integrated over in the definition 
of the PDF (for a specific illustration, see the introduc- 
tory sections of Ref. Q). In spite of this, the standard 
treatment is usually sufficient if one is concerned only 
with inclusive quantities. It is especially important that a 
well-defined operator definition exists for the PDF, since 
the universality of this PDF is part of what gives pQCD 
calculations their predictive power. 

It is now well-known that, in some situations, the 
transverse momentum of the struck parton becomes 
important. Examples include the treatment of the 
transverse momentum distribution in the Drell-Yan pro- 
cess [61 , and treatments of scattering at very high energies 
(small Bjorken-a:) where there is no transverse momen- 
tum ordering of emitted gluons ^] . This has lead to the 
use of transverse momentum dependent (TMD) PDFs 
which have explicit dependence on transverse parton mo- 
mentum. These are also called "unintegrated PDFs" 
(UPDFs) because the integration over transverse com- 
ponents of momentum is left undone inside the definition 
of the PDF. There is still an integration over the minus 
component of momentum, so the factorization formula 
still involves approximations that violate over-all four- 
momentum conservation. However, this is not a serious 
problem for many applications. 

An important observation is that generalizing factor- 
ization to include TMD PDFs is not a simple matter of 
leaving the integrals over transverse momentum undone 
in the standard definition of the PDF. Indeed, providing 
consistent operator definitions for the UPDFs involves a 
number of complications. One problem is that the most 
obvious candidate definitions suffer from divergences that 



arise from gluons with infinite rapidity in the outgoing 
hght-like direction. These "hght-cone" divergences re- 
main even after an infrared cutoff is included (see ]^ for 
a review of this and related issues). A consistent def- 
inition, therefore, requires an explicit rapidity cutoff in 
the definition of the UPDF [3, [lO|. Furthermore, it has 
been shown that factorization with TMD PDFs is vio- 
lated for some processes [ll|, [I^l ■ For an up-to-date re- 
view of the different types of PDFs, see Refs. [I,[i3|- See 
also Refs. [11 [H [H, IB, M, [11 for recent work related 
to defining the TMD PDF. 

The standard collinear factorization formalism and 
formalisms that use TMD PDFs accurately describe 
a wide range of phenomena. However, some stud- 
ies [20, [2l|, [22, [23] have illustrated that, at least for calcu- 
lations which require detailed knowledge of final states, 
an exact treatment of over-all kinematics is needed. In 
such an approach, the PDF should carry information 
about all components of parton momentum, including 
both the transverse and minus components. The defini- 
tion of the PDF in this formalism should not involve an 
integral over any component of momentum. A PDF that 
depends on all components of parton four-momentum is 
therefore called "fully unintegrated" (or "doubly uninte- 
grated"). In this fully unintegrated treatment, the kine- 
matics of initial and final states are kept exact through- 
out the calculation. However, as discussed above, the 
standard kinematic approximations are needed in the 
derivation of the usual collinear factorization formalism. 
Hence, as with a treatment of TMD PDFs, the fully un- 
integrated treatment requires a new approach to factor- 
ization. Among other issues, a precise definition for the 
fully unintegrated PDF is needed. 

A formal derivation of a fully unintegrated approach 
to factorization was initiated in Refs. [5J, [2J]. In [2J|, a 
complete, fully unintegrated treatment of factorization 
was given in the context of an MCEG. To avoid com- 
plications with gauge invariance, the authors restricted 
consideration to a scalar-0'^ theory in six space-time di- 
mensions. There, the derivation of factorization utilized a 
nested subtraction scheme in which double counting sub- 
tractions were required to obtain the correct hard scat- 
tering coefficient at each order of perturbation theory. 
This is a generalization of the approach to factorization 
discussed in [2^ and is reminiscent of the Bogoliubov 
approach to renormalization. Performing the subtrac- 
tions consistently required the introduction in [2J| of a 
well-defined mapping from exact to approximate parton 
momentum. The derivation in |24| was for arbitrarily 
many jets in the final state. Though the proof was for 
a scalar theory only, the procedure for mapping exact to 
approximate parton momentum is likely to extend to the 
gauge theory case. 

In [5| , a complete fully unintegrated derivation of fac- 
torization for a gauge theory was given for the case of 
a single out-going jet. The focus in Ref. ^ was on the 
leading order (LO) contribution to the hard scattering 
coefficient (though the single gluon vertex correction was 



also calculated). Since the relevant Ward identities do 
not yet have explicit enough proofs for the non-Abelian 
case, the derivation in [3] only applies rigorously to an 
Abelian gauge theory. The resulting structure is, how- 
ever, highly suggestive of a similar factorization formula 
for the non-Abelian case (QCD). 

An important difference between the standard 
collinear formalism and the fully unintegrated formalism 
is the number of non-perturbative functions that appear 
in a given factorization formula. In inclusive deep inelas- 
tic scattering (DIS), for example, a total cross section 
calculation that uses the standard collinear factorization 
approach only involves a single hard scattering coefficient 
and a single PDF. Schematically, the structure of the 
standard collinear factorization formula is: 



J2C<»f^/p■ 



(1) 



Here, C is a hard scattering coefficient and fi/p is the 
usual integrated PDF for scattering from a parton of type 
i inside the proton. The symbol (g) denotes the usual 
convolution integral. 

By contrast, the fully unintegrated approach requires 
not only a fully unintegrated PDF (FUPDF) in the ini- 
tial state, but also fully unintegrated jet factors and soft 
factors to describe final state parton four-momentum. 
Thus, the fully unintegrated factorization formula found 
in Ref. [5| for scattering from a target quark in the pro- 
ton includes not one but three non-perturbative factors: 
a fully unintegrated quark PDF, a fully unintegrated jet 
factor, and a soft factor. The integrals that would nor- 
mally allow the final state non-perturbative factors to 
be simplified are left undone in the fully unintegrated 
treatment. The structure of the factorization formula in 
Ref. [5| for DIS is written schematically: 



C ig) Fg/p (g) J (g) S*, 



(2) 



where F, J, and S represent, respectively, the fully un- 
integrated quark PDF, the fully unintegrated jet factor, 
and the fully unintegrated soft factor. Collectively, we 
refer to the fully unintegrated non-perturbative objects 
that appear in the fully unintegrated treatment as "par- 
ton correlation functions" (PCFs). In Ref. [| the PCFs 
in Eq. ^ are given precise definitions as expectation 
values of field operators. Clearly, the proliferation of 
non-perturbative factors makes the fully unintegrated ap- 
proach much more complicated than the usual formalism. 
The advantage, however, is that it provides a much more 
complete description of the distribution of final states. 

Much of the work done so far with PCFs and the fully 
unintegrated approach to factorization has been rather 
formal, with very limited phenomenological applications. 
Therefore, the aim of this paper is to initiate a closer con- 
nection to phenomenology by directly applying the meth- 
ods of Ref. [3, [2JI to a calculation of the hard scattering 
coefficient for deep inelastic wide-angle jet pair produc- 
tion from a target gluon inside a proton. The partonic 



subprocess is then, ^*g -^ qq (see Fig.[T]). This is the sim- 
plest possible calculation of a non-trivial next-to-leading 
order (NLO) hard scattering process in the fully uninte- 
grated approach to factorization. Since gluon-mediated 
processes dominate at small values of Bjorken-a;, the re- 
sulting hard scattering coefficient can be used to param- 
eterize the PCFs. Conversely, with a particular model 
(or parameterization) of the non-perturbative factors, the 
fully unintegrated hard scattering coefficient calculated 
in this paper can be used to make direct predictions that 
are consistent with factorization. Since the hard scat- 
tering coefficient is to be determined, as in [2J|, from a 
sequence of double-counting subtractions, the calculation 
relies directly on the results of @, UM ■ The result of this 
paper, therefore, is a prescription that allows for a sim- 
ple and direct implementation of the fully unintegrated 
approach in calculations of gluon mediated DIS. 

Although the work of Refs. [a, !24| is sufficient to pro- 
vide the steps for calculating the order-g^ hard scatter- 
ing coefficient, it is worth emphasizing that a complete 
derivation of factorization for the production of two jets 
requires a systematic consideration of arbitrarily many 
collinear and soft gluons. Additionally, the Ward iden- 
tity arguments in Ref. [5| need to be extended to the non- 
Abelian case. Since the purpose of this paper is to begin 
to apply the results of Refs. 0, [2J| to phenomenology, 
we leave these important but more formal issues to fu- 
ture work. Additionally, we set aside for now the issue of 
the evolution of the PCFs, initial/final state showering, 
or the relationship to other integrated or unintegrated 
gluon distribution functions. 

The paper is structured as follows: In Sect. [TTl we es- 
tablish conventions and set up notation. In Sect. UTTl we 
outline the basic steps in the calculation of a fully un- 
integrated hard scattering coefficient. In Sect. llVi we 
go through the steps for low-angle scattering and write 
down the LO factorization formula, in Sect.|V]we discuss 
scattering with large transverse momentum. In Sect. IVII 
we go through the steps for obtaining the subtraction 
terms that allow for a correct evaluation of the wide-angle 
(high transverse momentum) hard scattering coefficient. 
In Sect. IVIIl we write down the NLO factorization for- 
mula. In Sect. IVTIIl we summarize the results and make 
concluding remarks. 
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Following standard conventions, we work in a frame 
where the total transverse momentum of the photon- 
proton collision is zero, and express the incoming proton 
and virtual photon as 



P=\P^ 



2P+ 



'XP^ 



2xP+' 



(7) 
Since we insist on using exact kinematics throughout, x 
is not the usual Bjorken-x variable but rather 



2x 



Bj 



-4^x^ 



(8) 



1 + Vi 

with XBj being the usual Bjorken-x variable: 



XBj 



2P -q 



(9) 



Additionally, it will be convenient to work in the Breit 
frame where a:P+ = Qj\f2. The momenta of the final 
state jets are l\ and li^ and their masses are M\ and M^ 
respectively. The target gluon is labeled by k. Momen- 
tum labels are shown in Fig. [1] 

The contribution to the hadronic structure tensor from 
the graphs in Fig. [1] is given by the expression 



Wl^^{P,q) 



d% 



d% 



d'^k 



(\M{hMM'') 



flly 



4tt J (27r)4 J (27r)4 J {2^)^ 

X {2Tr)H'^^'> (k + q-h-h). (10) 



II. SET-UP AND NOTATION 

The process under consideration is described by the 
Feynman diagrams shown in Fig. [TJ The bubbles repre- 
sent target and final state subgraphs. All calculations are 
done in Feynman gauge. In a more precise treatment, we 
would also need to include a final state soft factor, but 
this will not be critical at the order of perturbation theory 
that we consider here. (If we continue to higher orders, a 
soft factor needs to be included, as in Ref. 0.) We will 
usually express four-vectors in terms of light-cone vari- 
ables. If T^ = {V~^, V~,Yt) is a four-vector, then we use 



We use the letter M to denote amplitudes for the full pro- 
cess: 7*p — > 2 Jets -I- X, in contrast to A which will be 
reserved to denote amplitudes for partonic subprocesses. 
The subscript g on the hadronic tensor symbolizes that, 
in the present treatment, we only consider graphs with 
a gluon in the target. In our analysis, it will be conve- 
nient to treat the final state jet momenta as independent. 
Therefore, we leave the integral over the momentum- 
conserving (5-function in Eq. ([TO)) undone until the very 
end. Throughout all calculations in this paper we use the 
convention that a factor of e^/47r is pulled out front as in 
Eq. (fTO|) . The subscript j labels quark flavors, and any 
sum over flavors is assumed to be implicit. A structure 





(c) id) 

FIG. 1: Graphs contributing to gluon induced production of two jets at order-gj . 



function is obtained from Eq. (jlOp by making appropri- 
ate projections onto the electromagnetic vertices. For 
example, for the Fi structure function we have 



F,iP,q)^P^,W^''iP,q), 



where 



P„ 



-9t^u 



Q^Pf.P. 
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M2Q2 



(11) 



(12) 



III. STRATEGY 

To calculate the hard scattering coefficients in Fig. [T] 
in a way that is consistent with the methods of p, [2J| we 
take the following steps: 

1. For each graph, categorize the leading 
regions of momentum space. 

The leading regions in momentum space for 
a given Feynman graph can be categorized ac- 
cording to the pinch-singularity surfaces (PSSs) 
of the Libby and Sterman analysis [26|. The 
two leading regions contributing in Fig. [T] are: 
the target collinear region, i?i, where the struck 
quark (ki — k ~ I2 in graph (a)) is collinear to 
the target proton; and the hard (wide-angle) 



region, R2, where the outgoing final state jets 
are at wide angles relative to each another with 
large transverse momenta. There is also a leading 
region, i?i, where the target-collinear parton is a 
struck anti-quark with momentum k'l = li ~ k (see 
graph (d)). We will always assume that the gluon 
momentum, fc, is collinear to the target proton: 



fc- 



A2 



(13) 



Corrections for gluon momenta away from Eq. ([T 
are to be dealt with in higher order calculations. 
(For an elementary discussion of categorizing lead- 
ing regions in terms of PSSs see, e.g., Ref. [27|.) 

2. For each leading region, apply approxima- 
tions on lines entering the hard scattering 
subgraph, appropriate to that particular 
leading region. 

In the scalar theory this only involves replac- 
ing exact kinematic variables by approximate ones 
inside the hard scattering subgraph. Here it is im- 
portant to note that initial and final state momenta 
are not altered at all in the approximation. 

A well-defined mapping from exact to approximate 
variables is required to ensure consistency with fac- 
torization. Such a mapping for arbitrarily many 



outgoing jets was given in Ref. [2J] for a scalar the- 
ory. In a gauge theory, one also needs to project ap- 
propriate external parton polarizations, as is done 
for the external quarks in Ref. [5]. In this paper, we 
will need analogous projections to extract appropri- 
ate gluon polarizations. Using the notation of Q, 
the complete set of replacements in the hard scat- 
tering subgraph, including polarization projections, 
will be symbolized by the action of an "approxima- 
tor" on the unapproximated expression, Eq. (jlOp . 
Since a different set of approximations is appro- 
priate for each region, then each region is asso- 
ciated with its own approximator. For example, 
Ti symbolizes the approximator appropriate to re- 
gion Ri, and TiW'^'''{P, q) is a good approximation 
to Eq. (Uni) in region Ri. So, TlW^"'{P,q) is to 
be read as an instruction to replace the exact mo- 
menta in the integrand of Eq. (jlOp with approxi- 
mate momenta appropriate for region i?i, and to 
make appropriate projections on parton polariza- 
tions. Analogous approximators, Ti and T2, are 
associated with regions ^1 and i?2- (Each of these 
approximations/approximators will be defined ex- 
plicitly in later sections.) In each case, the ap- 
proximator makes substitutions that are good ap- 
proximations in the corresponding leading region. 
This notation will make the logic of the subtrac- 
tions straightforward. It is important for this logic 
that, although each approximation is only good in 
a certain region of momentum space, each is well- 
defined for all momenta. 

3. In each leading region, use Ward identities 
to disentangle soft and collinear gluons into 
factors that are identifiable as contributions 
to PCFs. 

In the calculation of this paper, a Ward identity 
argument will allow the gluon to be disentangled 
from the out-going quark when the struck quark or 
anti-quark is collinear to the target proton (regions 
i?i and ^1). In the resulting factorized structure, 
one of the factors will be identified as the zeroth 
order contribution to the fully unintegrated hard 
scattering coefficient. The other factor will be 
identified as an order-g^ contribution to a pertur- 
bative expansion of the fully unintegrated quark 
or anti-quark PDF. The factor corresponding to 
a contribution to the fully unintegrated PDF can 
then be seen to be consistent with the operator 
definition for the fully unintegrated quark PDF 
given in Ref. [5']. 

For a detailed discussion of the PCFs and their def- 
initions, see section V of Ref. [a|. Here we very 
briefly review the definition given there for the fully 
unintegrated quark PDF. 

In Ref. [g], the fully unintegrated quark PDF was 



defined in coordinate space to be 

= {P\^{w)V^ins)In^.,^^o\VoinsMO)\P)R- (14) 

The momentum space distribution is then found by 
Fourier transforming. Let us briefly describe each 
of the elements of Eq. P^ : i/j is the quark field 
operator, and \P) is the target proton state. Wilson 
line operators, V,V\l, are needed in Eq. pi)) to 
ensure exact gauge invariance. Hence, we define 
Vu,{ns) as the path-ordered exponential of gauge 
fields from coordinate u; to cxd: 

VUn) = 

Pexp(-igs dXn-A{w + Xn)\ , (15) 

where A is the gauge field and P is a path-ordering 
operator. The Wilson line direction is nearly at rest 
in the center of mass system 



(- 



p2/s 



%0i) ys--0. 



(16) 



This direction is chosen to give the most univer- 
sal derivation of factorization |28j . Furthermore, 
the fact that the Wilson line direction is nonlight- 
like ensures that light-cone divergences associated 
with gluons moving with infinite negative rapidity 
are regulated [3, [l3] ■ In an extension to the non- 
Abelian gauge theory, the exponent in Eq. (fT5|) will 
also need to include the generator for the gauge 
group in the fundamental representation. 

For exact gauge invariance, Eq. pi|) also needs a 
Wilson line segment, In^-wfi, at light-cone infinity 
to connect the other two lines [23|. The definition 
of this gauge link at infinity is. 



In 



w.O 



Pexp -igs / dz^A^(z) 



(17) 



where C is a contour that connects the points at in- 
finity. This operator yields a factor of unity for our 
calculations since we work in Feynman gauge. The 
combination of Wilson line operators in Eq. (fT4|) 
forms a continuous Wilson line that connects point 
to point w. 

In Eq. (fT4|) . fj, is the renormalization scale associ- 
ated with the usual ultra-violet divergences. The 
subscript R on the right side of Eq. P^ signi- 
fies that renormalized operators are used. Thus, 
the fully unintegrated quark PDF depends on the 
quark momentum, k, the proton rapidity j/p, the ra- 
pidity parameter yg used to parameterize the Wil- 
son line direction, and the renormalization scale fi. 
For a detailed discussion of these issues and the mo- 
tivation for Eq. (fT4|) as the definition of the full 
unintegrated quark PDF, see section VC of Ref. 



4. Perform all necessary double counting sub- 
tractions. 



The complete faetorized expression for Eq. pTI)l is 
obtained via subtractions as in [2J, [25| . We briefly 
summarize the steps here. 

In a first approximation, contributions from high 
transverse momentum jets can be neglected be- 
cause they are suppressed by powers of g1 relative 
to the LO contribution while for low transverse mo- 
mentum we cannot rely on asymptotic freedom. If 
this level of accuracy were sufficient, then we would 
only need to analyze regions Ri and ^i , which were 
already considered in detail in 5] . That is, we could 
use the approximator notation from step 2 to write 

M^7(P, q) = TiI¥7(P, q) + TiI¥7(P, q) + 

+o((^)Vr(^''?)i)+^(.9?), (18) 

where a > and A is a typical hadronic mass scale. 
We symbolize the contributions from regions -Ri 
and Pi by 

W^:f^^^{P,q)^T,W^^{P,q)^ (19) 

wi;::,^,{p,q)^f,w^^{p,q). (20) 



We know from step 3 above that these contributions 
to the hadronic tensor will ultimately be expressed 
in terms of order-g^ contributions to the fully un- 
integrated quark (anti-quark) PDFs. 

To include the contribution from large transverse 
momentum scattering, we first recognize that the 
remainder terms in Eq. (|18p are given exactly by 
the subtracted expression, 



remainder = (1 - Ti - Ti)I^f"'(P, q) 



(21) 



Up to order-g^ corrections, this remainder term is 
dominated by wide-angle jets originating from the 
produced quark and antiquark, i.e., from the large 
transverse momentum region, P2, in the graphs of 
Fig. [TJ In approximator notation, T2 implements 
approximations that are good in this region. Thus, 
we apply T2 to Eq. ((2T|l to get a good approxima- 
tion to the contribution from region R2: 



K'a^miP^ 9) =T2{l-T,- T,)Wg^-'{P, q). (22) 



The subscript, 7*5 — + qq, indicates the relevant 
hard partonic subprocess for region i?2- Instead of 
Eq. pS]) we then have, 



Wg^'^iP, q) = (Ti + Ti)W/g'"'(P, g) + (1 - Ti - T^)W^''{P, q) 

= (Ti + T^)W^%P, q) + 72(1 - Ti - fi)W^''{P, q)+0 






(P, q) + Wi^'^^^^iP, q) + I^,^.";^,,(P, q)+0 



\W^'^{P,q)\^+0{gt) 
^^\w^'^{P,q)\]+0{gt). (23) 



A demonstration that the errors in Eq. 
actually A/Q-suppressed appears in Ref. 
scalar field theory). 



corresponding hard scattering coefficient will be of 



The first term in the last line of Eq. (|23p is the con- 
tribution to the structure tensor from the reaction, 
7*p -^ Jq + X, i.e., the production of a single jet 
due to the hadronization of a single knocked-out 
quark. The second term is the analogous contri- 
bution for the case where the jet originates from 
a knocked-out anti-quark. Both of these terms are 
examples of the zeroth order contributions to DIS 
that were already treated in Ref. [3]. 

The third term in Eq. ([^5]) is the contribution to 
the hadronic structure tensor from the reaction, 



7 P 



Jn 



h 



X, i.e., from the production of 
two wide-angle jets that result from the hadroniza- 
tion of a knocked-out quark and anti-quark. The 



orde 



r 9 s 



The last two terms represent the sup- 



pressed remainder terms. 



The first three terms in the last line of Eq. (P5)) give 
the contribution from the graphs in Fig.[T]to gluon- 
induced DIS, up to order-g^ and power suppressed 
corrections. The first two terms will factorize into 
a zeroth order hard scattering coefficient, a fully 
unintegrated quark PDF, and a fully unintegrated 
jet factor, just as in Ref. Q. The subtractions in 
Eq. (|22p will then allow the third term to be writ- 
ten in an analogous faetorized form. The aim of 
this paper is to find explicit expressions for each of 
the first three terms in Eq. ((23|) in terms of fully 
unintegrated PDFs and fully unintegrated jet fac- 
tors. 



In the rest of this paper, we will go through the de- 
tails of the steps enumerated above, ultimately arriving 
at an explicit expression for the factorization formula in 
Eq. 



REGION Rn TARGET COLLINEAR 
REGION 



IV. 



A. Figure [TJa) 

We begin by considering the target coUinear region 
Ri, starting with Fig. [IJa). The steps followed in this 
section are essentially the same as those described in 
Ref. 0, but illustrated in detail for the specific case of 
the graphs shown in Fig. [T] Since the quark is target 
collinear, we will apply approximations that allow the 
unapproximated graphs to be factored into a zeroth or- 
der hard scattering coefficient, and a factor identifiable 
as an order-g^ contribution to an expansion of the fully 
unintegrated quark PDF. We restrict the detailed dis- 
cussion here to the case of a struck quark because the 
treatment of the anti-quark case follows analogous steps. 

The components of the struck quark momentum in re- 
gion Ri are of order 



fci 






Recall that 



ki = k — I2 



(24) 



(25) 



We define approximate momentum variables appropriate 
for Ri with a "hat," 



These are exactly the values of the parton momentum 
obtained in the parton model. Now we set up a sys- 
tematic symbolic/diagrammatic notation to describe the 
approximator, Ti. For region _Ri we should treat the elec- 
tromagnetic vertex as the basic partonic interaction — 
it should be regarded as a separate unit which (after ap- 
proximation) can be analyzed independently of the rest of 
the graph. Therefore, we graphically denote the target- 
coUinear approximation by inserting a solid red circle 
around the electromagnetic vertex as shown in Fig. [21 
The circle should be read as an instruction to replace 
momenta passing inside it by approximate variables ap- 
propriate for Ri. Additionally, at each intersection of a 
propagator line with the circle, there is a projection ma- 
trix needed to project out the relevant Dirac components 
for unpolarized scattering. We call the lowest order hard 
scattering subgraphs on the left and right sides of graph 
(a) Ti-i^Q ^{hjki). (They are rather trivial at this order 
since they only involve the electromagnetic vertex factor 
7^.) The appropriate target-collinear approximation (as 
determined in fSI) is to make the replacement 



n^ 



LO 



,aih,ki 



'PTKoALki)r^ 



n^oAh,h)^rjn^o,aih,h)ri. 



(27) 
(28) 



fci-(-g+,0,Ot 



h^iO,q-,Ot). 



(26) 



The placement of the projection operators, Vt and Vj, 
is represented in Fig. [2^ by the points where the parton 
lines intersect the red circles. The projection operators 
are the same as in [5| and can be conveniently expressed 
in terms of the approximate parton momenta: 

Explicitly, the unapproximated squared amplitude aver- 
aged over color for graph (a) is 



(\Maih,l2,k)fy'^ ^g^.TnTi 



Yj{h)r 



I'^J^h)!" 



*.p(fc,n, 



(30) 



where the color factor is Tr = 1/2 for an ^[/(S) gauge theory and to is the quark mass. The factors $„p(fc,P) and 
J(l\) and J{l-2) represent, respectively, the incoming target subgraph and the final state jet subgraphs. These are 
represented by the bubbles in Fig. [51 The red circles in Fig.[21are to be read as an instruction to make the replacements 



7^ -^ V.n^V.j, 7" -^ VtI^Vt, 



(31) 



after which the amplitude becomes 

(2Q2)2 

-2\ M'' 



ri(|A/„(/i,/2,fc)r 



Tr 



Q' 



/i7^lil Tr \kj{li) 



A7*9-«(fi,fci) {TTlr-Jih)]}\g',TRTr 



Tr 

7+ 



/i 



rJihh" 



l^J{l2)Y 



1 



^i — m 
1 



$Kp(fc,P) 



(32) 





FIG. 2; Figure [ija) with the apphcation of the target coUinear approximation for region Ri. The graph on the right-hand side 
of the arrow shows the separation into the factors of Eq. (I32p . 



where 



A'''i^'}{h/ki) 



-Tr 



"iirk 



7 ill 




(33) 



is just the 7*9 — > q amphtude squared and averaged over spin. (RecaU that in our convention a factor of e^ has aheady 
been factored out front in Eq. PH)) .) In Eq. ([5^ we have used the approximator notation discussed in Sect. IIIIl Ti 
acts on the amphtude in graph (a) by replacing it with the approximate expression on the right side, appropriate for 
region i?i. 

The trace in the second factor in Eq. ([5^ comes from expanding J' in a basis of the Dirac algebra 

J{1) ^Js + 7^ Jm + ^''"^M- + 75 J5 + 7^75 Jm5, (34) 

and recognizing that the dominant term in unpolarized scattering is J^^ {li)j~^ ■ This is the term projected out by Vt 
and V.J in Eqs. ()31|) . The minus component of J7{li) is determined from the trace 



J-(l,)^-TT[^-Jih)]. 



(35) 



In the second line of Eq. (|32p . we have used braces to help distinguish the basic elements of the factorized expression. 
Once the integral over I2 is included, the last factor in braces is identifiable as an order-^^ contribution to the fully 
unintegrated quark PDF, Eq. HH). The simple partonic vertex Eq. ([55)1 has been pulled out in front. We have 
extracted an over-all factor of 1/Q^ to give the conventional normalization of the 7*9 — > q amplitude. The factorized 
form in Eq. ([5^ resulting from the action of Ti on Eq. (jTU]) is illustrated by the graph on the right side of the arrow 
in Fig. [21 To simplify notation, the factored bubbles in Fig. [2] implicitly include the traces in Eq. ([32]) . Hence, to 
simplify later expressions we define 



Jih)=TT[^-J{h)] 



J{l2)^Tr[^-J{h)] 




(36) 
(37) 



B. Figure [T](b) 



For gauge theories, the situation is more complicated than, for example, the six-dimensional scalar-;/)'^ theory used 
in Rcf. ^245 because, in addition to graphs like Fig. [Ija), we can also have leading contributions from graphs like 
Fig. [Ijb) where there is both a target coUinear quark and a target coUinear gluon attaching to the hard part (see 





FIG. 3: Figure [TJb) with the red circles showing the apphcation of the target colhnear approximation. The graph on the right 
hand side of the arrow shows the separation into the factors of Eq. (I42|) . The gluon is shown attaching to an eikonal hne 
represented by double lines. 



the left side of the arrow in Fig. [3]). The extra coUinear gluon can be disentangled with the help of a Ward identity 
argument as we now illustrate explicitly. The unapproximated expression is 



{\Mt,{hM,kfY'' ^ glTn^y 



I'Jih)!^ 



YJih)!" 



^.p{k,P). 



(38) 



Following step 2 of Sect. IIIIl we draw solid red circles around the hard parts as shown in Fig. [3] On the right side 
of the cut, we make the same approximation as in graph (a). However, on the left side, we must first dot k into the 
hard bubble and multiply by n'^/ins ■ k — ie), as dictated by the prescription in section XB of Ref. Q. Inside the hard 
circles, the prescription of Ref. [5| is to replace momentum variables by 



h = iO,q-,Ot), 



' = '^'°'°*'' 



-k^ 



q+li 



k+ -U 



,0,0, 



(39) 



Note that ki — k — I2, which is consistent with the requirement that approximate momentum is conserved. The solid 
red circles therefore tell us to make the replacements 



7 



1 



ti-^- 



r 



■ k 



-rj% 



1 



/ 1 — k 



l^Vj, 



7 



VtI^Vt- 



(40) 



Note that it is necessary to use the precise definition of the replacements in Eqs. (|39l40p to remain consistent with 
the factorization formalism of Ref. |5] . The contribution from graph (b) to the complete amplitude thus becomes 



Mh 



fiiy 



(2Q2)2 



Tr 



fli% 



h-% 



rh 



Tr 



kj{h) 



k — i> 



-Tr 



iiJ{h)Y 



^.p{k,P). (41) 



A Ward identity argument now allows the coUinear gluon to be factored away from the hard part; substituting 
% = -{li- %) + /i into Eq. (gT]) leads to 



1 



(l^bl')' 



(2Q2)2 
■2\ 1^'^ 



Tr 



7''/i7^^il Tr %J{h) 



Ar'<!^i{h,ki) 



{TT[^'Jih)]}i~g^,Tji 



■k 



■k- 



-Tr 



-Tr 



/l^(^2)'- 



-^^(^2)7'' 



$Kp(fc,P) 

^.p{k,P)\. (42) 



C. Figures [T](c) and [T](d) 

The remaining graphs in Fig. [3] follow exactly similar steps as graphs (a) and (b) . For graph (c) we have 
Ti(|Af,|^)'"' = 

-Tr ^{^, )7"J(/2 



A'^'i^'i{hM) {Tr[7-J(/i)]}^.g,^Tj 



■k 



ris ■ K + le 



4 \ <c^ — m 



$«,(fc,P) , (43) 
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and for graph (d) 



Ti 



AT*9-«(/i,fci) 



■2\ A"' 



{Tr[7-J(/i)]} g^r^ 



n; 



Us ■ k — ieUs ■ k- 



-Tr 



ze 



'-J(h) 



^.pik,P)\. (44) 



D. The sum of graphs (a)-(d) 

Adding graphs (a-d) results in a factorized structure which we can express in terms of diagrams as 



T^{\M{h,h,k)\'y'^ ^-1^ 




(45) 



The graphs in the first two factors in Eq. (|45p are the same as in Eqs. (|33l35p . They do not depend on l2- Hence, we 
can push the I2 integral in Eq. (jlOp through to integrate over the last factor in parentheses and obtain the single real 
anti-quark contribution to the fully unintegrated quark PDF, 



^i/?(^i,^) 



(2^) 



' 9'sTr<Ti- 



1 



Us ■ k 



-Tr 



^Jihh" 



— ni 
1 



nP 



Us ■ k 



1 
Tr 



- m 
nP 



fZ Ti a ' tZ 



Tr 



+ 



'-Jih] 



^^p{k,P). (46) 



The superscript (Iq) symbolizes that these diagrams represent the contribution to the fully unintegrated quark PDF 
due to the hadronization of a single real anti-quark in the final state. Note the double lines in the last factor of Eq. ([IS]) 
denoting the eikonal propagators of Eq. (j46p . Equation (|46p is precisely what arises from an expansion of the operator 
definition in Eq. (|14p for the fully unintegrated quark PDF. (See Ref. [ll, la] for more details on writing down Feynman 
rules for the quark PDF.) In Eq. (j46p and from here on out, we drop the explicit ie terms in the eikonal propagators 
to simplify the expressions. Given models (or fits to data) for J7 and $, we can explicitly calculate the contribution 
from Figs. [1] to the fully unintegrated quark PDF in perturbation theory using Eq. ()46p . Of course, perturbative 
methods are only valid for a large enough kf. Near the core of region i?i (fci with almost exactly vanishing minus 
and transverse components), the fully unintegrated quark PDF should be modeled or fit to data. Additionally, there 
is freedom in how one chooses the precise value of j/s so long as y^ « 0. 



E. Region Ri: Target anti-quark 

If we consider the case that the target coUinear parton is an anti-quark, then instead of Eq. (|24p we use kinematic 
approximations appropriate to region ^1: 



where (see graph (d) of Fig. [J) recall that 



fc;~(Q,AVg,A) 
k[ ^h~k. 



(47) 



(48) 
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We define the approximate variables analogously to Eq. ((26)) for the target quark case. This time it is the anti- 
quark which gains a large minus component after being struck by the virtual photon. So, using an inverted "hat" to 
symbolize approximate variables for region i?i, we have 

k[ = i-q+,0,Ot), l2^i0,q-,0t). (49) 

The remaining steps are exactly analogous to those of the treatment of the quark target in the previous subsection. 

F. Summary: Factorization formula for 7*p ^ 1 Jet + X scattering. 



We can summarize the results of this section by writing the explicit factorization formula for the contribution to 



the hadronic structure tensor from 7*p - 
from a knocked-out quark (anti-quark). 
(5-function, we have 



Jg(g) + X scattering. (The subscript q{q) on J indicates that the jet arises 
Using Eq. (|15)l inside Eq. (jlOp and using the ^i integral to evaluate the 






■J,+x{P,<l) = ± 



d^ki 1 
(2^Q2 



ylT*9-9(?i,fci) 



■2\ A"' 



J(;i)F,/p(fci,P). 



(50) 



We have changed variables so that the integration in Eq. (|50p is over struck quark momentum, fci. Furthcrnrore, we 
have dropped the (Ig) subscript on the fully unintegrated PDF in Eq. 1(501) because, in general, we need to include 
all graphs that contribute to the quark PDF in addition to those shown in Eq. ((46|) . Note that the non-perturbative 
functions, J{li) and Fg/p{ki,P) are evaluated using exact momentum. The procedure for evaluating Eq. 150]) can be 
summarized as follows: 

• Obtain parameterizations of J(/i) and Fg/p(ki, P) from models or fits. 

• Evaluate the LO hard scattering coefficient, A'' '^^*, using Eqs. ((^^ to obtain li and ki. 



The analogous expression for i?i is 



I'p- 



'Jn+X 



{P.q) 



3 

47r 



d^k[ 1 

(2^Q2 



Ut* 



9-9(^2, fci)| J{h)Fq/p{k[,P). 



(51) 



Given a particular model or parameterization of the fully unintegrated jet factor, and the fully unintegrated quark 
PDF, Eqs. (|50l5ip allow a direct calculation of the first two terms in Eq. ((^5() . Together, they give the complete 
contribution from the graphs in Fig.[l]to LO hard scattering. Away from the core of regions Ri, one can also estimate 
Fq/p by directly calculating Eq. (|46p . given knowledge of <i>(A:, P). For low- angle scattering we should also include the 
order-gg virtual correction to the hard vertex. This was already worked out in Ref. 51 so we do not repeat the steps 
here. 



V. REGION R2: LARGE TRANSVERSE 
MOMENTUM 

The next step is to consider region R2 , where the virtu- 
ality of ki is large and Eq. (|26p is a poor approximation. 
However, k still has nearly target-collinear momentum 
given by Eq. P^ . and ^1 and I2 are highly boosted rel- 
ative to one another. Therefore, we should regard the 
7*(7 — > qq subgraph as the basic partonic subprocess; af- 
ter making approximations appropriate for R2 , it should 
be possible to factor the elementary on-shell amplitude 
for 7*5 — > qq scattering. Therefore, the approximated 
graph in region i?2 will be symbolized diagrammatically 
by a dashed blue rectangular box enclosing the 7*3 -^ qq 



subgraph as shown in Fig. (H The box should be read as 
an instruction to replace all momenta inside it by the ap- 
propriate approximated momenta for R2 (to be discussed 
below). Furthermore, at each intersection of a parton line 
with the dashed-blue box, there needs to be an appro- 
priate projection onto Lorentz/Dirac components. The 
purpose of this section is to establish exactly what these 
approximations must be such that the resulting expres- 
sion is consistent with the factorization formalism that 
was set up in Refs. [^, l2J| . 

Returning now to Eq. ((ID)) , we decompose the square 
modulus of the amplitude into the contraction of sub- 
graphs 



J 



{\M{hM,k)\- 



= Tr 



Jih)H^L^ih,l2,k)J{l2)H''ji'' ih,l2,k) 'S>^,.'{k,P) 



(52) 
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(a) 



(b) 



FIG. 4: Dashed blue boxes symbolizing the application of the approximation appropriate for region R2 to Figs, [ija) and (b). 
Analogous graphs are needed for Figs.[T][c) and (d), though we do not show them here explicitly. 



The Ti-L/R are now the hard subgraphs for region i?2j i-C, 
the subgraphs enclosed by dashed blue boxes of Fig. S) 
To simplify notation, Dirac indices will be left implicit. 



A. Approximations on kinematic variables 

Inside the dashed blue boxes, we replace the exact 
parton momentum variables with approximate momen- 
tum variables appropriate for the production of two 
wide-angle jets with large transverse momentum. We 
use the approximation scheme proposed in [2J] applied 
here to the special case of two out-going jets. The 
choice of a good approximation scheme is probably not 
unique, though once one is chosen it should be main- 
tained throughout the rest of the calculation. The ap- 
proximate momentum variables for region R2 will be de- 
noted by a tilde placed over the corresponding unapprox- 
imated variable. 

Enforcing partonic values for the approximate mo- 
menta amounts to requiring 



momenta that satisfies these constraints is \2A 



fc=(fc+,0,Ot), 



;? = o. 



li = 0. 



(53) 



The true values of kj, fc~, l\ — M^ and Z| = M2 there- 
fore parameterize the deviation of the approximate kine- 
matics from exact kinematics. In addition, we require 
that both the exact and approximate momenta obey to- 
tal four-momentum conservation: 



k + q = li+l2j 
k + q = li+l2- 



(54) 
(55) 



~k-i q+ +':'+'/ AoA 

\ 2/r 212 J 


(56) 


fltl - - \ 

\2l^ J 


(57) 


h- ~,l2~h.t , 


(58) 


li^, =li,, - k,/2, \2,t = h,t - kt/2; 


(59) 


[f =Zf - k-12, f^ =l2 - k~/2. 


(60) 



where 



In the limit of small kj and small k , we have \j^t 



1 



3,t 



and I; ~ L . Again, the choice of approximation scheme 
in Eqs. (|56ll60p is not uniquely determined from Eqs. ((551 ^ 
I55p . Other equally good approximation schemes rnay 
have advantages, but the one used here and in Ref. [2J| 
is sufficient for our purposes. With the constraints of 
Eq. (|54l55p we also have li_t = —h,t- Given the exact 
values of initial and final state momentum, we can write 
the inverse transformations 



l,.t + 



j-t -l'3,t 
P 

1+ _S:t 



2' 



^J=n 



2L 



E^.^ 



Here j — 1,2 label the out-going quark and anti-quark 
jets. We therefore have a one-to-one mapping between 
exact and approximate momentum variables for region 
i?2. It will also be useful to define approximate versions 
of the partonic Mandelstam variables: 



~s = {q + kf, i={l2-kf 



u={li-kf. (61) 



A choice of a mapping between exact and approximate 



Our main goal is to apply a sequence of approximations 
that allow the unpolarized, on-shell square modulus of 
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the partonic amplitude to be factored out of the complete 
sum of graphs in Fig. [TJ Part of the approximation is, of 
course, to replace h^h, and k by their approximate values 
li, I2, and k, defined above. That is, we will ultimately 
need to make the replacement 



where 



K%ih,l2,k) 



nl%{h,h,k). 



(62) 

However, in a gauge theory we also need to project out 
the relevant Lorentz/Dirac components of the external 
lines, in a way that is appropriate for unpolarized scat- 
tering with massless on-shell partons. The steps for these 
projections will be discussed in the next few subsections. 



B. Target gluon polarization 

To factorize the fully unintegrated gluon PDF, we will 
first need to project out the relevant polarizations of the 
target gluon. It will be convenient to work at the level of 
the amplitude, decomposed into upper and lower blocks 
as shown in Fig. [5l The sum of graphs contributing above 
the dotted line will be called U and the sum of graphs 
below the dotted line will be called C. The steps for 
projecting gluon polarizations will follow closely those in 
Ref. [sol for the standard integrated PDF. Decomposing 
the amplitude in this way, we have 



M^'{h,h,k)=U^:{h,l2,k)g-^ C^'{k,P), 



(63) 



Here we have explicitly extracted the numerator, g'^'^, 
of the gluon propagator. We will neglect the electro- 
magnetic indices and the momentum arguments in this 
expression when convenient. 

The projections on target gluon polarizations should 
allow the lower part of the graph to be separated into a 
factor that is consistent with a definition for the fully un- 
integrated gluon PDF. To this end, we note that, in the 
center of mass frame, the target bubble is highly boosted 
in the plus direction so the components of the lower bub- 
ble have relative sizes given by, 






o 



A 



(64) 



Since none of the lines in the upper part of the graph are 
collinear, all of the components of U are comparable: 



u- -^W 



'Ut. 



(65) 



Given Eqs. (|64l65p it is tempting to keep only the term 
U~ C'^ in (p5|) . On a graph-by-graph basis, however, this 
leads to the appearance of super-leading contributions 
(i.e., contributions to structure functions that vary as 
Q^ - see, e.g., Ref. [3l|)- These super- leading contribu- 
tions correspond to scalar and longitudinal polarizations 
whose total contribution cancels exactly in the sum over 
all graphs, at least for the case of a single target gluon. 
Therefore, they should be systematically removed graph- 
by-graph. We do this, following the Grammer and Yen- 
nie treatment treatment in QED [32| , by first writing the 
gluon propagator as a sum of two terms which we call 
"iiT-terms" and "G-terms," 






kt'n': 



fj,l^ 



k'^n': 



(67) 
(68) 



We use ns in these definitions because this suggests a fully 
unintegrated gluon PDF density that is closely analogous 
to the quark density defined in Ref. [5J. The iiT-term in 
Eq. ((66)) vanishes exactly in the sum over all contributions 
to the upper bubble in Eq. (j63p . so we drop it in our 
calculations. Therefore, we replace Eq. ([63]) with 



U-{h,l2,k)C^ik,P). (69) 

In the last equality, we define £ as a short-hand for 
the contraction of C,^' with a G-term. Only the trans- 
verse components are unsuppressed in Eq. ()69p . Recall- 
ing Eqs. ()13l64p we can verify this directly by checking 
the relative sizes of each combination of indices 






(70) 



(71) 



Thus, after the replacement in Eq. ([55]) . the plus and 
minus components for a given graph are suppressed rela- 
tive to the transverse components. (Recall that ys ~ 0.) 
Note that it is important that the contribution from the 
K-ierva vanishes exactly in the sum over all graphs in or- 
der to avoid uncontrolled errors, so we must use k rather 
than k in Eqs. ([68l[67|. 

As part of the wide-angle approximation, we drop the 
power suppressed terms in Eqs. (|70l7ip and make the 
replacement 

M(;i,Z2,fc)-^^z^^(;i,Z2,fc)>c,(fc,p)- 
j 

Y,Y.^\hM,k){tt,,y {tt,,y c,{k,p). (72) 

The sums over i,j are only over transverse momentum 
components. In Eq. (j72p we have introduced the trans- 
verse polarization vectors 

(etj)i = (0,1,0,0), (etjf^ (0,0,1,0). (73) 

Using Eq. (|72p in the squared amplitude, summed over 
final states, we have 



r.iJ-l' — ft-A"' 



A'' 



_Q^.v 



(66) 
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U 



FIG. 5: Basic amplitude separated into an upper part U and a lower part C, see Eq. (|63|l 



\M{hM.kf^ E E (E"'(^i'^2,A:)Z^^"'t(;i,Z2,fc)] (e,,,r(e,,,r(e,,,,r'(eM'r'(E^»(^'^)4(fc,^))- 

i,j,i'.j' s,s' \Xj / \Xt / 

(74) 
Here ^j^ and J^x ^^^ sums over the final states for the jets and the target respectively. We can use the fact 

that, in an unpolarized cross section or structure function, W{li, I2, k)W ^(^i, I2, k) is diagonal in j and j' to re-write 
Eq. dUD as 



{\M{WMM^Y^ ^ (iEE"-''(^~i'^~2,fc)z^-'tai,f2,fc)) iY,T.^''^G^'''' c^{k,p)ci,{k,p) 



Xj j 



Xt j' 



^EETr [jih)nr{hM,k)J{h)n'^' {hM.k)] {e,,,ne,,,yFl%{k,P). (75) 



To symbolize the gluon projections described above, we use a projection matrix T'?-, , defined to implement the 



(76) 



replacement in Eq. ((72| . That is, the replacement in Eq. ((72)) is expressed as 

M{h,h,k) ^ Wihj2,k)r^^,a'{k, P). 



C. The fully unintegrated gluon PDF 

The square modulus of the amphtude in Fig.[5]is factored into two parts in Eq. (jTS]) . The first factor is simply the on- 
shell squared matrix element for scattering off a transversely polarized gluon, averaged over transverse polarizations. 
The second factor 



^S(^'^) = EE^''"^'''''^-(^'^)4'(^'^) = E^''"^''"'*--'(^'^) 



(77) 



Xt j' 



j' 



is what we identify with the lowest order correction to the unpolarized fully unintegrated gluon PDF. This suggests 
the following operator definition for the fully unintegrated gluon PDF: 



Fg/p{k,P) 



V^ / dw+dw d^Wt _,fc.^ 



j 



{2ny 



e-*'=-'"(p| [us^^gf^^iw)] [n,^,g-'Ho)] \P). 



(78) 



The sum over indices j involves only the transverse components, g is the gauge field strength tensor. The factorized 
structure of Eq. (|75|1 is shown graphically in Fig. [6] where the graph beneath the horizontal dotted line now corresponds 
to Eq. ([75]) . and the part above the dotted line corresponds to the first factor in parentheses in Eq. ([75)1 . The cross 
drawn at the end of the target gluon is meant to symbolize that it is only the G-term that is retained in the target 
gluon propagator in Eq. ([77]) . In a non-Abelian gauge theory, Eq. (I78p will need to be modified by the insertion of 
a Wilson line operator, similar to the definition in Ref. jj|, but these modifications won't effect the LO structure in 
Fig. El 
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FIG. 6: The square modulus of the amphtude after keeping G-terms and projecting transverse components of the gluon 
propagator. The disconnected upper and lower parts represent the factors in Eq. (|75p . The cross on the gluon in the lower 
bubble represents the contraction of the G-term with the lower bubble in Eq. 



D. Final state jets 



Next we make projections onto relevant Dirac components for the outgoing jet lines, appropriate for region R2. The 
final state jet bubbles can be expanded, as in Eq. (j34|) . in a basis of the Dirac algebra. The dominant components 
in these expansions are most easily found by working in coordinates where the +2 axis lies along the direction of 12- 
Then I2 has a large plus component and li has a large minus component. The dominant contributions for each jet 
bubble are then 






-7— —ir/ 2* 11/^(^2)* 1*2- 

{2sY 



(79) 

(80) 



Since we restrict to unpolarized scattering, we have dropped polarization dependent terms. In Eqs. (|79l80p we have 
used the approximate momenta to write the substitution in a frame independent way. Now it is clear what the Dirac 
projection operators are for region R2: 



V~, 



UiU 



V, 



tP^- 



Part of the approximation for i?2 is then to make the substitutions 
Equivalcntly, we may write 



l^J'\h) « ^Tr [l2J{h) 
l^J^'ih) - ^Tr yhjih) 



Making these substitutions in the squared amplitude, we have 



\M{hM,k)i 



^Tr 
= Tr 



Vj^Jih)Vjn^L^^{h,l2,k)Vi^J{l2)VjH''ji'' ih,l2,k) 



Ii'Hl''^{li,l2,k)l2'H''ji'^ {li,l2,k) 



iTr 

2S 



hJih) 



Finally, for convenience we define the following notation 



1 



4(^i) = ^Tr[/2:7(/i) 



JiSl2)^Y.TT[hjil2) 




JiM 



(81) 
(82) 

(83) 
(84) 



<^..K'ik,P) 

^TT\hj{l2)\}'i>^,^'{k,P). (85) 
2s L J J 

(86) 
(87) 
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The subscript on J indicates which vector is used for projecting Dirac indices. Note the difference in normahzation 
from Eqs. (|55I57)) . 



E. Factorized expression for wide-angle 7*p -^ 2 Jets + X scattering 



We now have aU the steps needed to factor the on-sheU 7*3 -^ qq matrix element from the rest of the graphs in 
Fig. [1] Combining aU of the approximations discussed in this section, appropriate for region i?2, amounts to replacing 
the hard scattering subgraphs by, 

H^''=(/i, /2, k) -. ViTi'ii^ihM, k)Vj^V^^,, 



K^'-ih.h, k) -> vt,,vjv!i:'{lui2, k)vj^. 



(89) 



The approximator, T2, appropriate for region R2 acts on the squared amplitude by making the complete set of 
replacements from this section in the graphs of Fig. [T]and Eq. ([5^ : 

T^ {\M{hM,k)\'Y = Tr \Vi^J{l{)Vin^^\hM,^ViJ{h)Vin''ji\hMM] V^,-K'^''''{k,P) 
= ?EE(^*.j)'(^*.J')'Tr \hK'\hj2,k)h'H''/ihj2M Jj^{h)Jj^{h)Fg^p{k,P) 



s jj' 



2\ f^" 



Ar'9^m{h^l2,k) \{h)\{l2)Fg/j,{k,P). (90) 



In the second line we have repeated the steps that give Eqs. (|75l85p . and we have dropped the small components in 
Eqs. (|70l7ip . The ordinary unpolarized on-shell matrix element has been factored out of Eq. ([90]) . 

To summarize, Eq. (I90p is the factorized expression for the graphs in Fig. [T] that is a good approximation in the 
region of momentum space corresponding to two high transverse momentum jets. The matrix element for massless 
on-shell 7*17 — > qq scattering with transversely polarized gluons has been factored out of the rest of the graph. The 
last line in Eq. (|90p can be interpreted diagrammatically as 



T; 



(|Af(/i,/2,fc)|' 



^y 




■hM 



■'hCi. 




(91) 



To obtain the complete NLO result, we need to avoid double counting by subtracting from Eq. ((90|) the LO result of 
Sect. HVl This is the topic of the next section. 



VI. DOUBLE COUNTING SUBTRACTIONS 



Recalling Eq. pS)) . the NLO contribution to the hadronic tensor is 

W^^P^g, (^' 9) = ^2 W^r (P, q) - T2T1 W^7 (P, q) - T2f, W^" {P,q). (92) 

The first term in Eq. ((92)) is found by using the amplitude, Eq. (|90|) . in Eq. (fTO]) for the structure tensor. The second 



(93) 



term is the subtraction term corresponding to region Ri and is found by directly evaluating 

T2T,{\M(h,l2,k)\^y\ 
Similarly, the subtraction term for region i?i (target collinear anti-quark) is found by evaluating 

^2^l(|Af(^l,^2,fc)|')''^ 

With the careful formulation given so far of the approximations for each region, it is now straightforward to obtain 
the correct subtraction term by directly applying the approximators, Ti,Ti,T2, in Eqs. (|93l94p . In this section we 
go through the steps of applying these approximations and we derive explicit expressions for the subtraction terms. 
Returning again to Fig. [I] we consider each graph separately. 



(94) 
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(a) 



ib) 



FIG. 7: The sequence of approximations needed to obtain the subtraction terms corresponding to Figs. [Tla)-(b). Analogous 
approximations are also need for Figs. [Uc)-(d) though we do not show them here. 



A. Subtraction term for region Ri 

1. Figure\^a) 

We start by calculating Eq. ((93|) . beginning with the box diagram, Fig. [T^, and the corresponding amplitude 
Eq. ([50)1 . We first apply the approximator Ti as it is defined in Sect. If VI and then the approximator T2 as it is defined 
Sect|Vl As in Sect. IfVI the application of the first approximator is represented by soUd red circles. The application 
of the T2 approximator is represented by a dashed blue box enclosing the entire ^*g -^ qq block including the solid 
red circles (see Fig. [T]). Inside the dashed blue box, parton momenta are replaced by the approximate values in 
Eqs. (|56l57l58p (each denoted by a tilde). Additionally, there is a projection at the intersection of each parton line 
with the dashed box. Inside the solid red circle, parton momenta are replaced by the approximate values in Eq. (j26p 
(each denoted by a hat). Again, there are projections at each intersection of a parton line with the solid red circle. 

The parton lines on either side of the electromagnetic vertices intersect the solid red circle so we make the replace- 
ments 



as prescribed by Ti. Next, the external legs of the 7*5 
the further replacements prescribed by T2, 



7^ ^ Vjj^'Vj, 7" -> VtI^Vt, (95) 

qq subgraph pass through the dashed blue lines, so we make 



Vj-f^Vj 



r 



1 



fi 



7 



m 



VtI^Vt 






(96) 

(97) 



There is a projection for each intersection of a parton line with the dashed blue boxes. Additionally, the momentum 
fci has been replaced by fci (Eq. ([5^)1 ) inside the boxes and the masses of the internal partons are set to zero. Making 
these substitutions in the square modulus of the amplitude leads to 



T2T, {\Ma{h,l2, fc)|') = 5? TnTr ^j-'VTVj^J{h)ViVjrVj i^^j -I'Vj^JihWi^l' [j^ Vt 

(98) 
Repeating the steps of the last two sections allows Eq. (j98| to be written in a more compact form with clearly 
recognizable factors: 

T2T^{\Ma{hM.k)i 



2(2S)2 



EE(^*.^)'(^M)^Tr 



fVThVj^^^Vj ( ^ ) 7^/: 



(27 



fc'^- 



{Xr [hJih]] } {Tr [hjih]] } Fg/p{k, P) 



= 2^^ 



''il-/^k]^a{ll,l2,k)\—TT\hjill) 



2S 



2S 



:Tr 



hj{l2)\\Fg^p{k,P), (99) 
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where in the second hne we define the factor 



^aihyh^k) 



(2Q2)2 



EE(^*.^-)'(^M)^Tr hh 



Tr 



^,3 



^■'^)^''-'(^ 



(100) 



Evaluating the contribution from Fig. [3)Ja) to Eq. (^0)) and then subtracting Eq. ([M)l yields the contribution from 
Fig. [Ha) to the first two terms in Eq. 



(\Ma{ll,l2,k)\^y'^ - T2T1 [\Ma{ll,h,k)\' 



aX 



S-»99 



Q-iM.k) 



■2\ A"^ 



^aV■\^ hi k) 



A-r*?-9(/i,fci) 



■2\ t^^ 



JlSh)\{h)F,,^{k,P) (101) 



One can verify, by repeating the steps above for Ti , that Fig. [Ija) does not contribute to the third term in Eq. (|92p . 
Thus, the factor in braces in Eq. (jlOip is precisely contribution to the NLO fully unintegrated hard scattering coefficient 
arising from the application of the approximations and subtractions to Fig. [Ija). Apart from Sa(/i,^2,^) (defined 
in Eq. (|100[) ) Eq. (|10ip only involves ordinary on-shell squared amplitudes for parton scattering. Graphically, the 
contribution to the hard scattering cocfhcient (in braces in Eq. ()101|) ) is 





(102) 



The second factor here is the subtraction term. The approximations that lead to Eq. (|10ip allow the hard vertex to be 
factored out of the rest of the subtraction term. Thus, we may write the graphical representation of the subtraction 
term in the simplified form: 




Sa(^l,Z2,fc) 




(103) 



In the subtraction term, after the vertex factor is extracted as in Eq. (jl03p . the singular behavior at fci -^ ki is 
contained in the factor Sci(^i, I2, k) which can be calculated directly by first writing it as 



2 rp 

Ea{h,l2,k)= ,'^^ g^pTr 



(2Q2)2 



kill 



Tr 



where we define the transverse tensor 



9L 



-9kp 



k^nj^p + nj,^,kp 
k ■ nj 



h { ^\ I'-hY ^ ' 



nj = (0,l,Ot). 



Evaluating Eq. (|104p yields 



=a(ti, f2, k) = - ~ (2fci • h) \Ali-k 



H 



k ■ nj 



k-nj 



(104) 



(105) 



(106) 



Note that 2^ depends on the exact momenta, li,l2, and k, via the dependence of the approximate momentum variables 
on exact momenta as defined in Eq. (|26I56II5S1) . 
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2. FiguresU}(b)-(d) 

We may also go through the steps for Fig. [Ijb). The sohd red circles and blue dashed boxes denoting the sequence 
of approximations are shown in Fig. [7]^b). Starting again with Eq. psp . we first apply the approximator for region 
i?i. The red circles tell us to make the replacements 



7 



/i-^ 



Vjj" 



1 



ii ~ k 



rvj, 



VtYVt- 



(107) 



Additionally, in accordance with our prescription for region Ri, we dot k into the left-hand circle and multiply by 
n^ / {us ■ k). So we make the further substitution 



P,77" 



h~% 



rvj 



n, • k 



Vj% 



h~i 



rvj 



(108) 



The outer rectangle tells us to apply the approximations for region i?2 to everything inside it. Therefore, on the 
left-hand side we make the substitutions 



■k 



Vji 



1 \ 71^ -7 



li — k 



l^Vj -^ V'Vi 



n 



■k 



1 



li — k 



I'^VjVi, = 0. 



(109) 



The left-hand hard subgraph in graph (b) exactly vanishes because 7^? , involves only transverse components (see 



3] 



Eq. [72]) whereas Ug has no transverse components. A similar result holds for graphs (c) and (d). Thus, only Eq. (|10ip 
contributes to the region i?i subtraction term. 



B. Subtraction term for region _Ri 



The steps for calculating the anti-quark subtraction term, Eq. ([M|) . follow exactly analogous steps to what was 
already done in the previous subsection for region i?i, so we do not repeat the details. Instead of Fig. [ija), it is 
Fig. md) that gives the only non- vanishing subtraction term in region Ri . The analogue of Eq. poip is 



(^\Md{liM,k)\' 



fJ.V 



ToTi 



{\Md{ll,l2M- 



fJ.1^ 



A}^-^^{k,h,k) 



2X A"' 



Ed{hj2,k)l Ar'^^'i{l^,k[) 



.\{h)\{h)Fg,j,{k,P). (110) 



The factor in the subtraction term analogous to Eq. (|106p is 



- r/ / i.^ alTn -, J - ~ {2h ■ k',)(2h ■ nj) {2h ■ h){2k[ ■ nj) 



k- nj 



k ■ nj 



(111) 



VII. COMPLETE NLO CONTRIBUTION 

Gathering the results of the previous section together, we may write down the order-g^ contribution to Eq. pop : 



W^!f -iP a) 



dHi 



47r J (27r)4 J {2ttY J {2ttY 



d^k ~ 



W!;.!fg^^^ih,h, k) 4(/i) .\ (h) Fg/p{k, P) {2n)H^^\k + q-h- h). (112) 



This is the desired factorized form for the third term in the last line of Eq. 
coefficient is, from Eqs. (|101lllOp . 



The fully unintegrated hard scattering 



W^l^g^^.ihM^k) 



A^'a^i^ihMM 



■2\ A"' 



Ea{h,l2,k) A'/*9-9(/i,fci) 



■2\ M!' 



Ed{hj2,k){ U7*9-9(;2,fci) 



/ji> 



(113) 



20 



In terms of diagrams, the hard scattering coefRcient is 

'1 

'U r- 





^aih,h,k) 




^d{h-,h,k) 




(114) 



The steps for evaluating the integral in Eq. (|112p are as foUows: 

• Obtain values of ^i, I2 and k directly from parameterizations or models of Jj (Zi), Jj {I2), and Fg/p(k, P). 



Determine the correct approximate values of parton momentum {Ii,ki,l2,k[,li,l2,k) using Eqs. (|26l49l56ll58p . 

Evaluate Eq. pi3p using these approximate momenta. The factors, Eaih^h^k) and 'E,dih,l2,k) are to be 
directly evaluated using Eqs. (|106lllip . The squared amplitudes in Eq. (|114p are just the usual ones obtained 
using Feynman graphs, but evaluated with the appropriate approximate parton momentum. Calculating these 
amplitudes requires that we make a projection on electromagnetic indices to obtain a particular structure 
function. For completeness, we give explicit expressions for the squared partonic amplitudes, with projections 
on electromagnetic indices, in App. |^ 



The first term in Eq. (|113p contains singularities at fci —> fci and k[ — > k[ which are exactly canceled by the 
subtraction terms where the singularities are contained in the factors, Sa(/i,/2,fc) and S(i(Zi, ^2, ^)- This is similar 
to what occurs in the usual subtraction approach to factorization. However, in the standard integrated approach 
to factorization, the subtractions involve generalized functions (e.g., ^-functions and "+" -distributions), so that the 
cancellation only makes sense after an integration over final states. By contrast, the fully unintegrated hard scattering 
coefficient Eq. (J113p is just an ordinary function that can be evaluated for any final state momentum. The subtractions 
of Ref. [a, |2J| that ensure consistent factorization occur point-by-point in momentum space. Using Eq. (JA1IIA6P in 
App. [XI we have verified numerically that the cancellation o f singularities takes place. 



VIII. SUMMARY AND CONCLUSION 

With the results of this paper we are able to directly 
calculate the gluon-induced contribution to the unpolar- 
ized structure functions to order g^ in a fully uninte- 
grated approach. The key results are the LO contribu- 
tions, Eqs. (|50l5ip . and the NLO contribution Eq. ITT^ 
to the hadronic structure tensor, Eq. (|10p . The LO and 
NLO results should be summed to accurately describe the 
overlap of the different regions of momentum space. Since 
we have used the subtractive method of Refs. [2J, l25j . 
then we can be confident that the errors are A/Q sup- 
pressed point-by-point in momentum space. 

Evaluating Eqs. (|50l51llOp only requires the expres- 
sions for the matrix elements obtained from ordinary 
Feynman graph methods (see Eqs. (|Aip - (|A6p ). and the 
explicit expressions given for the factors So(Zi, I2, k) and 
^d{hj2,k) in Eqs. (|106lllip . Finally, the approxi- 
mate variables used to evaluate the different factors (jAlf 
IA6ll06lllip are given by the mappings in Eqs. (|26l39l56f 
1611) which relate exact to approximate momentum vari- 
ables. 

Since there is a large overlap of regions i?i, Ri with 



R2, the subtraction terms need to be dealt with very 
carefully. As an example of how things could go wrong, 
we could imagine that we approximate k in the numer- 
ator of the subtraction factors, Eqs. (|106lllip . by k 
from Eq. (j39p . This could simplify the expressions and, 
naively, it might seem to make sense given that the pur- 
pose of the subtraction terms is to remove double count- 
ing from i?i and ^1 where fc is a good approximation to 
k. However, the subtraction terms are used in regions 
of momentum space that are far from the core of the Ri 
and Ri regions. Thus, replacing fc by fc in Eqs. (|106lllip 
yields large modifications to Eq. (|114p in certain regions 
of phase space. Alternatively, one could imagine try- 
ing to change the definition of k in Eq. ([5^ . However, 
the definition of k that we have used is rigidly fixed by 
the requirements of the factorization formula in Ref. [5| , 
with a well-defined fully unintegrated quark PDF. We 
can thus see how the lowest order factorization formula 
and the corresponding definitions for the fully uninte- 
grated PDFs places important restrictions on the higher 
order calculations. 

We have so far considered only the 7*3 — > qq par- 
tonic subprocess. To complete the NLO treatment re- 
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quires that we include the somewhat more comphcated 
7*9 — > gq partonic subprocess, so the result of this 
paper is not quite a complete NLO treatment. How- 
ever, this paper provides a concrete illustration of the 
techniques developed in Ref. [^, [2J] applied to order- 
g^ scattering, and provides a starting point for con- 
sidering more complicated diagrams. Furthermore, the 
graphs in Fig. [1] dominate processes at low- a; and are 
the only graphs included in some MCEGs. Thus, the 
results presented here are already of phenomenological 
interest. The treatment in this paper can be imple- 
mented directly in Monte Carlo calculations similar to 
what was done in Ref. [23J, |33| where it was shown that 
fully unintegrated PDFs are needed to obtain an accu- 
rate description of the distribution in final states. Using 
the results of this paper, one has direct access to the 
fully unintegrated gluon PDF. Furthermore, the tech- 
niques described here are likely to be compatible with 
existing methods ([13, [SH, [H, His, 36, 37]) in the de- 
scription of initial and final states; a description of the 
non-perturbative bubbles, $ and J', can likely be ob- 
tained by utilizing these other approaches. 

In the large- a; limit, the standard kinematic approxi- 
mations of collinear factorization can affect even totally 
inclusive quantities, and a fully unintegrated formalism 
may be needed. The set-up in this paper naturally takes 
large-x kinematic effects into account. An alternative, 
simpler approach was developed in Ref. [3a | where cross 
sections where written in terms of the usual integrated 
PDFs of the standard collinear factorization theorems, 



but such that large-x effects from initial and final state 
masses were taken into account. It is possible that this 
approach can be extended to the order-f^^ jet produc- 
tion calculation presented in this paper in the large-x 
limit, thus simplifying the evaluation of Eq. (|112p . (In 
the large- a; limit, however, it is not sufficient to keep only 
the diagrams of Fig. [U and quark-induced processes also 
need to be included.) 

The natural next step is to extend the calculation il- 
lustrated in this paper to the case of single gluon emis- 
sion. That is, for the partonic subprocess, we also need to 
consider diagrams for 7**7 — > gq scattering. The calcula- 
tion will follow steps very similar to those presented here. 
However, it will be somewhat complicated by the need to 
take into account subtractions that correspond to contri- 
butions to the soft factor and the jet factor. Therefore, 
the soft factor will need to be explicitly included every- 
where. 
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APPENDIX A: PARTONIC AMPLITUDES 



Structure functions like Eq. (jlip are calculated by making projections on electromagnetic indices using suitable 
combinations of g^'^ and P^^P^ . Here we write down the explicit expressions obtained by projecting with g^^'^ and 
pt^P'^ on the partonic squared amplitudes in Eq. (|113|) . For the LO partonic amplitude the result is 



9i^u 



P^.P. 



A-y'i^Hh^ki) 



■2\ M" 



v4''*9"*9(/i,/ci) 



■2\ A'l' 



= -2Q^ 



= 0. 



(Al) 

(A2) 



9^A\A'i''i-"'{h,K)\ = -2Q 



fll^ 



P^.P.{ U'/*9-?(/2,A:i) 



(A3) 

(A4) 



For the NLO partonic amplitude the result is 



ffpi^ 



^7*9-99(^1, ^2, A:) 



2\ A"' 



-4=9s Tb[- + j r^ 

,u t tu 



(A5) 



and, 
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P^^P. 



A'>'9'-<}']{l^j2,k) 



2\ A"' 



= -2g',TR{{2k-P) 



2/1 P 2l2- P 

i u 



t u , 



+ ^^ {2~sQ^Ml - (2P • fi)( (2P • fci)(2^; • h) + < (2P • fc^) ) - (2F • l2){ (2P • l[){2h ■ h) -u{2P- h) )) . (A6) 



We have verified numerically that Eq. (IA6p is non-singular at the core of regions Ri and Pi . The desired unpolarized 
structure functions are obtained by taking appropriate combinations of g^t^ and P^^Pi, as in Eq. (jlip. 
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